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1. Introduction 

One of the many interesting developments to arise out of Maldacena's conjecture []I|- 
^ concerns the study of large A^ QCD. While it is not known how to construct the dual 
gravity theory for SU{3) gauge theory with six quark flavors, nor even the dual for any 
pure SU{N) theory, it is hoped that the gravity duals that are known will share many of 
the properties of everyday QCD. The main properties we seek are asymptotic freedom in 
the ultraviolet and conflnement in the infrared. 

In four space-time dimensions, the Maldacena conjecture was originally applied to A/" = 
4 supersymmetric Yang-Mills [0 and then to other theories that have less supersymmetry 
but are still conformal [§|-14|. Since the /3-function is zero for these theories, they do not 



behave like ordinary QCD. 

Witten has proposed supersymmetric Yang-Mills at finite temperature as a model 
more suited for comparison to QCD |T^. A d dimensional euclidean gauge theory at finite 
temperature is equivalent to a theory with d — 1 noncompact directions and a Euclidean 
time compactified on a circle. At large distances the theory acts like a nonsupersymmetric 
euclidean Yang-Mills theory in d — 1 dimensions. Hence, one expects to find area law 
behavior and a gap. The gravity dual for finite temperature super Yang-Mills is a nonex- 
tremal black hole in AdS space. Using this, Witten was able to argue that the bulk dilaton 
wave equation has a discrete spectrum, implying a gap for the boundary theory. He also 
demonstrated that there is an area law for the space-like Wilson lines. It was later shown 
that the finite temperature theory exhibits magnetic screening |l^,l7|. There has also been 



a mini-industry comparing dual gravity results to lattice results |T^ , with some reasonable 
agreement [|T9|-p7|. 



However, in the UVthe finite temperature theories effectively become supersymmetric 
with one extra dimension. Thus, the coupling does not run, but remains a free parameter. 
Related to this problem is that the QCD scale is set by the temperature instead of by 
dimensional transmutation. So the glueball mass scale is roughly the temperature T, 
while the string tension, with the free parameter, is g'^NT'^. In a QCD theory where the 
scale is set by dimensional transmutation, one expects the glueball mass scale to be directly 
related to the string scale. 

We thus seek another gravity dual to better describe QCD. One such candidate is the 
type model proposed by Klebanov and Tseytlin [|2^ . Ironically, this model was considered 
a toy model when first developed ||2^,0. Following an original suggestion by Polyakov 
|31j] , Klebanov and Tseytlin argued that one could construct the dual of an SU{N) gauge 
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theory with 6 real adjoint scalars by stacking N electric D3 branes of the type model on 
top of each other. The type model has a closed string tachyon that lives in the bulk, but 
no open string tachyon that lives on the branes. The closed string tachyon couples to the 
five form field strength, which then drives the tachyon to a nonzero expectation value. 



In [^ it was shown using the equations of motion derived in p^ that the effective 



coupling between two heavy quarks has the desired logarithmic fall off in the UV 

Qym ~ ^j , (1-1) 



where (i is the energy scale. A possible IR solution was given in [^, but after a closer 
inspection, one finds that for generic tachyon potentials, the UV solution does not connect 
to this particular IR solution. In [^| a different IR solution was given which corresponds 



to an IR fixed point with infinite coupling and no confinement. For a given class of tachyon 



potentials, such a fixed point could conceivably connect to the t/F point in [p5 |. 

At first glance, the presence of this IR fixed point seems sensible since, as was shown 
in 



3^ , the two loop contribution to the /3- function is positive. Hence one might expect a 
Banks-Zaks fixed point ||3^. Presumably, quantum corrections can push the value of this 
fixed point to infinite coupling. 

On the other hand, one could argue that quantum corrections lift the masses of the 
adjoint scalars. Hence this theory should behave much like pure SU{N) gauge theory, with 
some extra massive states. Using this reasoning, one should expect the type theory to 
confine. 

In fact, using only some modest assumptions about the tachyon potentials, we will 
show that the classical dual gravity equations of motion have generic solutions in the IR 
which exhibit i) a linear quark- antiquark potential ii) magnetic screening and iii) a discrete 
glueball spectrum. The glueball mass scale is related to the string scale by m ~ ^ jN^I'^, 
with no dependence on a free parameter. If one assumes that the tachyon reaches a fixed 
value in the /i?, then there exists IR solutions where the curvature diverges, but where 
many, if not all, a' corrections fall off exponentially. 

In section 2 we review the type model and discuss a W solution with a running 
coupling constant. We assume here that the tachyon potentials satisfy a minimum number 
of properties, but are otherwise generic. In section 3 we discuss various IR solutions. We 
argue that in order to reach the IR solutions discussed previously in the literature requires 
plenty of fine tuning, and even then, it is not guaranteed that the UY solution of section 
2 can match to these solutions. We then discuss a generic class of IR solutions and argue 



that the UV solution is adjustable enough to match over the range of these IR solutions. 
In section 4 we examine the validity of these solutions beyond the classical limit, where in 
general we find that string corrections swamp the leading order term. However, we find 
some solutions where many a' corrections are exponentially suppressed. In the IR these 
solutions approach a conformal transformation of the product space i?-(- x M4 x 5*5, with 
the tachyon flowing to a constant value. This leads us to speculate that there is an exact 
solution to the full cr-model whose asymptotic behavior matches the asymptotic behavior of 
these classical solutions. In section 5 we discuss the physical implications of the generic IR 
solutions in section 3. Given one inequality, which is satisfied by the special solutions found 
in section 4, we argue that two heavy quarks have a linear potential at large distances. We 
further argue that there is magnetic screening and a discrete spectrum. We also calculate 
some glueball mass ratios using the WKB method. The results differ somewhat from the 
finite temperature results and the errors are larger due to the complexity of the equations. 
For the low lying masses we can make reasonable estimates for lower bounds of mass ratios 
which are consistent with lattice computations. In section 6 we present our conclusions. 



2. The Type Nonsupersymmetric Model. 



The type model |]2^j33] has a closed string tachyon, no fermions and a doubled set 
of R-R fields |]5^,0 , and thus a doubled set of D branes [^ . Because of this doubling of 
the R-R fields, one can relax the self dual constraint on the 5-form field strength. Hence 
one can have D3 branes that are electric instead of dyonic. The low energy world- volume 
action for A^ parallel electric D3 branes is SU{N) QCD with six real adjoint scalar fields, 
but no fermions. Hence, there is no supersymmetry and the coupling will run. There is no 



open string tachyon [Q, so there is no tachyon in this QCD model. 



Klebanov and Tseytlin argued that the massed squared of the closed string tachyon 
gets a positive shift from the background 5 form flux [^. With a background flux the 
tachyon potential is not symmetric under T -^ — T, so a large flux can drive the tachyon 
expectation value to a nonzero value that is basically independent of the tachyon bare mass 
term and its quartic coupling. The tachyon field is a source for the dilaton, thus unlike 
the A/" = 4 case, the dilaton expectation is not constant. 



One then makes the following ansatz for the metric |^^ 



ds'^ = e^^ (e^^-^'^dp^ + e-^^dxl + e^^'^^d^f) , (2.1) 



where (/>, ^ and r] are functions of p only. The equations of motion then reduce to a Toda 
like system with an action [E8i 



(2.2) 



S = j dp 1^2^i^2^iT2-5r)2-y(<^,e,r7,T) 

and a constraint 

\4>^ + ^^' + \t^ - 5f + ^l*/*, ?, ^) = 0. (2.3) 



That the left hand side of ( p.3| ) is a constant follows from the equations of motion. The 
fact that the constant is zero follows directly from the dilaton equations of motion and the 
Einstein equation. Q is the total D3 brane charge which is proportional to A^, T is the 
tachyon field and f{T) is a function given by [^ 



/(T) = 1 + T+-T2. (2.4) 

This describes the tachyon coupling to the five form fiux. The bare tachyon potential is 
the negative of fi'(T) 

g[T) = \t^ - XT' (2.5) 

where we have included the quartic piece of the potential. The coefficient A of the quartic 
coupling has yet to be determined unambiguously, so for the purposes of this paper, we 
will treat it as a parameter. Actually, the details of /(T) and g{T) are not that important, 
except that we require that f{T) have a minimum at some finite value of T and that g{T) 
is positive at that value. We will often consider the special case where g'(T) = when 
/'(T) = 0. 



The equations of motion derived from the action in (|2.2| ) are 



I + -^(T)e5'^+^^-^^ + 2 



Q 



2 



e-2^ = 



(2.6) 



rj + -g{T)e^^+^^-^'' + Se"^^ = 
^ + -^(T)e^'^+5«-5'? = 

For large Q, g{T) plays a secondary role to /(T) and so the tachyon expectation value 
is determined by setting f'{T) = 0. As a first approximation, we may assume that the 



tachyon is a constant T = Tq. If Tq = then the solution reduces to the Af = 4 solution. 
When T is nonzero, then all fields are coupled and there is no known analytic solution. 
However, we can attempt to find approximate solutions that are valid in the UV and IR 
regions. 

An asymptotic solution valid in the UV was given in [^ . As a first approximation 
one can take the dilaton field to be relatively constant, at least compared with ^ and 77. 
Assuming a constant cj), the equations for ^ and rj can be solved exactly, at least in the 
near horizon limit. In this case we find 

e^ = Cip e^ = Csp^/^ 

1 c^/^ 1 202 

2S(T„)^a.* + ^ -1 = p.7) 

59(r„)^e4* + 1^ - 6 = 0. 

One can easily check that this satisfies the constraint equation in ( |2.3| ). If we plug this 
back into the metric, we find that the solution is still AdS^ x 5*5, but the curvatures of 
the two spaces no longer match; 5*5 now has smaller curvature then AdS^. In this case the 
Ricci scalar for the total space is proportional to 

R ~ 9{To)e^^. (2.8) 

Using the ^ and rj solutions as inputs, we can go back and find an approximate solution 
for (J) in terms of p. Using the ansatz es*^ = Co{log{p/po))°', and plugging this into the 
equation of motion for (p in (|2.6|) , we find that the ansatz is a leading order solution if 
a = —1 and Cq = —4C2/{g{To)^/Cl). po is an integration constant and we assume that 
po >> 1 in order that the gauge theory length scale is much greater than the string scale. 
Setting p = u~^, and using the lowest order solutions for Ci and C2 from (|2.7|), we learn 
that the leading order behavior for the coupling is 

1 Q^^(To)(log^)' 

e-'^ = ^^ = ) °^ . (2.9) 

g^M 1024v/27(7i;y 

One can easily check that to leading order in 1/logu., the constraint equation is still 
satisfied. We can also estimate the range of validity for this solution. Computing the 
leading order corrections to Ci and C2, one finds that 



We can also compare the terms in the potential that depend on the tachyon. Since 

^^(To)e^^+^?-^- ~ ^^ Q'f-\To)e-^^ ~ ^, (2.11) 

our solution with constant T = Tq and /'(Tq) = is valid so long as \og{u/uo) >> 1. 
While we are treating A as a parameter, we should note that it is crucial that g{T) > 
when f'{T) = 0. Hence we assume that < A < 1/2. The lower bound on A is so that the 
complete tachyon potential is bounded below. 

The metric in the large u limit is 

9{To)\og^^\u^ 2Q \ log^J II 1^ log^; ' 

(2.12) 
Hence we can trust the classical dual gravity solution only if g{To) << 1, since 
log(u/wo) >> 1- However, for generic values of A, g{To) ~ 1 when /'(Tq) = 0, hence 
it would seem that the classical result is not particularly trustworthy. As it happens, the 
situation is not completely hopeless [0. While the curvature for ( |2.12| ) diverges in the 
large u limit, the Weyl tensor actually falls off as 

\C\ ~ -^. (2.13) 

logw 

This is a consequence of the conformal invariance in the large u limit, where the space-time 
approaches that of AdS^ x 5*5, with divergent curvature. Because of this fall off in |C|, the 

o 

a' correction is of the same order as the classical term. Hence, while quantitative results 
might not be exact, the classical dual gravity solutions could capture the true qualitative 
behavior for this SU{N) gauge theory. 

It would seem that the dual gravity calculation failed its first test: the prediction in 
( p^) is that e~'^ has a log squared dependence instead of the linear log behavior found 
in perturbative Yang-Mills. However, the physical coupling is determined by finding the 
potential between two heavy quarks. Using the Wilson line computation of |^,0, one 
finds that the quark potential is given by [p5[ 



V ^ - ^^f/" L«Lo (2.14) 

where Lq is some length that can be adjusted to be much longer than the string scale. 
Hence, the effective coupling between a heavy quark and its antiquark does appear with 

6 



the expected log dependence. In deriving (|2.14|) we have substituted R 



16 



g(To)log(u/uQ) 



into the quark potentials found in in [^|36| . A perturbative Yang-Mills computation gives 
a quark potential proportional to g'^N. Since the perturbative Yang-Mills coupling behaves 
like g'^N ~ l/log(Lo/L), we find that the dual gravity computation results in an effective 
coupling with the desired log dependence. The actual coefficient is related to the one loop 
/3-function. Unfortunately, it is clear that the coefficient is model dependent and we don't 
know enough to actually compute this using the gravity dual. Even if we did, we should 
not necessarily expect to get the correct answer since stringy corrections will be of order 
1. 

As we move away from the ?7F point at w = oo, the tachyon will begin moving away 



from Tq. Solving the equations of motion to the next leading order we find that ^2 
a'{To) 1 , ^/log(-logp) 



B 

+ 




logp 

(2.15) 



logp 

1 , /. \ 1 ^ /^log(-logp) 

77 = - log(4p) - ^-_ + O ( ,^ 2 ^^ 



2 logp V log P 

The constant B in the (j) expansion is an integration constant that can be removed under a 
conformal rescaling. Notice that Q only appears as an overall factor in front of e^. We also 
notice that the sign of g'iTo) determines whether or not the tachyon expectation value 
is driven toward or away from zero as one moves away from p = 0. There is also the 
possibility that g'iTo) = 0, in which case the tachyon expectation value is unchanged as 
the system varies. 

Another interesting point concerns the ratio of the two and one loop contributions to 



the /3- function [32]. To next to leading order the effective coupling ea*^ is given by 



1 



e'^ ^ 77tY2 • (2-16) 

^OgU - (g + g(Vo)2 ) loglogW 



Let us compare this to perturbative Yang-Mills, where the coupling, up to two loop order 

1 

log U — 2^- log log U 



IS, 

1 



gvM^ ~ 1;——^ (2.17) 



and where 61 and 62 are the one and two loop contributions to the /9-functiona. For an 
SU{N) gauge theory with six adjoint scalars, the ratio is 62/26^ = 3/16. Hence we find that 
the sign of 62/26^ coming from dual gravity is model independent, and that the minimum 
value of the ratio, and the one that comes closest to the perturbative result, is 7/8 when 
9'{To) = 0. 



3. Connecting to the Infrared 

While it is satisfying that a reasonable UV solution exists, it is not immediately clear 
what the IR behavior is like. Numerical simulation is very difficult with four coupled 
nonlinear differential equations, so while useful, it does not immediately lead one to the 
answer. 

To better explore the situation let us consider a greatly simplified model where 
/'(To) = g{To) = g'{To) = 0. In this case the tachyon is constant and the other three 
fields decouple from each other. Obviously the solution for the dilaton is ^ = aop + (f)Q. 
The other equations of motion reduce to 

^ mf (3.1) 

rj + 8e~^^ = 0. 

The AdS^ x 5*5 solution is ^ = log{y^2 / f {Tq)Q p and f] = ^ log(4p). However, there can be 
other solutions to ( p.l| ). Consider the case where the solution asymptotically approaches 
the AdS^ X 5*5 solution in the limit p ^ 0. Each second order differential equation should 
have two integration constants. Two of these are fixed by choosing the singularity at p = 0. 
These generalized solutions are given by 



e = log(v/277(7i;yQp) + J^a^p^ 

(3.2) 



ry = -log(4p) + ^hnp 



2n 
•J-nP 

n=l 

2n 
WnP 

n=l 



where ai and 61 are free parameters and the higher a„ and 6„ are determined by recursion 
relations. 



A more careful calculation of the effective coupling from the Wilson loop shows that there 
are no additional log log « terms. We use the convention that 61 is negative and 62 is positive. 



We can also consider corrections to the IR (p -^ oo) AdS^ x 5*5 solution. In this case 
the solutions have the form 



e = logi^/yJiT^Qp) + J2^np- 



n=l 



1 ^ (3.3) 

77 = -log(4p) + }_^dnP ", 



2 

n=l 



where ci and (ii are free parameters and the higher Cn and dn are computed via recursion 
relations. It is straightforward to show that a UV solution with generic values for ai and 



6i does not connect to the IR solutions in (|3.3D for any value of ci and di. To see this 



note that ^ and r] both satisfy an invariance equation, 

U' - 7^e-2^ = Ci - 5r)2 + 20e-4'' = C^. (3.4) 

^ Jy-i-o) 



A quick calculation shows that for the [/F solutions, Ci = 2ai and C2 = — IO61. However, 
the IR solution has Ci = C2 = for any ci and di. Hence, unless ai = bi = 0, the 
UV solution does not connect to this IR solution. To see what solution it does connect 
to, note that if ^ >> 1, then the nonlinear term in the equation of motion is small, 
likewise for the equation for 77. Hence for large ^, the leading order behavior for ^ and 77 is 
^ = aip + Co + 0{e~'^"^P) and 77 = a2p + 'no + C>{e~'^"^P), with ai, 02 > 0. Note that these 
solutions are quite generic since they involve 4 integration constants. The constants of the 
motion are Ci = ai/2 and C2 = —ba^. Hence we see that the t/^F solutions connect to 
the IR solutions with ai = 2y/ai and a2 = \p2h\. If we include the a^ term from (/>, then 
the constraint condition in (^]^) implies 

From the C/F point of view, after fixing the coupling and the singularities at p = 0, we are 
left with three integration constants. The constraint reduces this to two. Adjusting these 
integration constants adjusts the values of a^, i = 0,1,2, which also satisfy a constraint 
relation. Notice that one possible transformation is the conformal rescaling, where 

i{p) ^ ^M - log(t) ^(p) ^ ^itp) - \ log(t) 4>{p) ^ 4>{tp). (3.6) 

Obviously, under this transformation the a^ rescale to ct^ — > toi, and so the ct^ can be set 
arbitrarily close to zero, the limiting result being the usual IR fixed point of AdS^ x 5*5. 
Aside from this conformal transformation, we are left with one independent transformation 
to the integration constants. 
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Let us turn to the more difficult case wliere /'(Tq) = g'{To) = 0, but (/(Tq) 7^ 0. 
From tfie last line of ( |2.6| ) we see that T is constant, hence the problem is reduced to 3 
coupled nonlinear differential equations. As was argued in the last section [P3| , P^ , there 
is a [/F solution whose leading behavior is given by (|2.15|) . Adjusting ^, r] and (p to all 



have singularities at p = takes care of three integration constants. Another integration 
constant is varied under the conformal rescalings in ( |3.6| ). This leaves two integration 
constants. These are of a similar type as the previous example, although they are entangled 
due to the coupling of the equations. The leading correction for ^ is 5ap^ while the leading 
correction to r] is ap^. These also generate corrections to (p. There is another integration 
constant appearing at order p^ which can be adjusted to satisfy the constraint in ( |2.3|) , 
leaving only one adjustable integration constant. 

Notice that since the t/Feffective coupling is (7^ ~ l/(logM), these subleading terms 
look like instanton contributions. In the perturbative regime, we should expect instanton 
contributions to be powers of u~^^'^ = p'^^'^. Hence in the large A^ limit, no such terms 
should appeara. Since all fields lie in the adjoint representation, it is possible that there 
are fractional instantons. However, these usually require fermion zero modes which are 
conspicuously absent here. 

Another source for u~^ corrections are scalar masses. These could be bare masses or 
masses generated by quantum corrections. The ability to adjust the integration constants 
on the type gravity side could then correspond to the freedom to adjust the bare masses 
on the field theory sidea. The powers of u could also just be an artifact of the strong 
coupling expansion. In any event, since this theory is not conformally invariant, there 
does not seem to be any a priori reason why they should be set to zero. Hence, we will 
assume that these coefficients are fixed, but generic. 



Let us now try to match to IR solutions. One such solution was discussed in [|25 
Consider the ansatz where e~^^ is small compared to the other terms in V. This corre- 
sponds to a small curvature for the 5 sphere. Dropping this term, one can now find an 
exact solution to the equations of motion that satisfies the constraint. The solution is 

e^ ^ dp'/^ e^ = C2p^/^ e« = /^ p, (3.7) 

with the relation 

10(2/(To))i/4C|-9(7(To)v^(7o = 0. (3.8) 



Unless of course A'' = 3! 
I thank I. Klebanov for a comment suggesting this. 
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Comparing all terms in V, one has e^'^~^2^~^'^ ~ p~2^ ^-2S, ^ ^-2^ -^yi^ e~^^ ~ p~^'^'^. 
Hence this solution is valid for large p. From ( p.7| ), the coupling blows up as p ^ 00 and 
after substituting p = 1/u^ the metric is 

C2 remains as a leftover integration constant. As in the [/F, the curvature in the IR 
is small if 5'(To) << 1. However, there is no reason to expect (/(Tq) not to be of order 
unity. Hence, as in the [/K, the string corrections are of the same order as the classical 
contribution. 

While e^ blows up for this particular IR solution, one can easily see that the effective 
coupling found from the heavy quark potential reaches a finite value. Defining a new 
variable v such that 

1 (^2f{% 



9 y 3Q 



1/2 

^/^/^ (3.10) 



the metric in (|3.9| ) is then 



5r(To) \ ^;^ II QH^ 



rf«^ = ^ PV + v^dxl + ^ ^^ V2 ' ^ ^^5 • (3.11) 



3 (To) 



Prom this metric, we see that R^ = -pf-r, and so the heavy quark potential is ^5 



720 7r3 

Hence, we see that this particular solution is a conformal fixed point in the IR. 

A possible cause of this behavior is that the SU{N) gauge group has been Higgsed to 
U(^l^ri-i j-jy ^YiQ adjoint scalars. Naively, one expects a repulsive force between the branes 
because of the extra R-R field. In the large N limit, the branes could be pushed apart, but 
still maintain the 5*0 (6) symmetry. The effect would be to have the branes smeared out 
over some region with spherical symmetry and finite size. In the UV, well above the Higgs 
scale, the gauge group is unbroken and the coupling runs. In the infrared, the Wilson loop 
will probe down into this smeared region and see the effect of the broken gauge symmetry, 
hence the quark potential behaves coulombically. 

However, this ignores the role of the tachyon, which one might expect to lead to 
attraction between the branes. From the perturbative side this seems reasonable, since 
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quantum fluctuations would give masses to the scalar fields, removing the flat directions 
away from the unbroken gauge theory. There is a Coleman- Weinberg potential that is 
unstable |^ , but this requires fine tuning to get rid of the scalar mass. 



Connecting this IR solution to the desired UV solution is problematic. Just as in the 
uncoupled case, the generic IR behavior for p -^ oo has the form 

(p ^ aop + (pQ 

e ~ «iP + ^0 (3.13) 

T) ^ 02/0 + rjo. 

The constraints on q:^ are 

a^ > 0, -al + -al - hal = 0, (3.14) 

where the latter constraint also insures that 5ai > ^ai + |q;2, thus insuring that all terms 
in the potential V have an exponential falloff. The UV solution is connected to the IR 
solution in (^]^) by varying integration constants such that the cti in (|3.13| ) are set to zero. 



Prom the IR side of things, we see that the at can all be set to zero by an infinite rescaling. 
But we do not want to do this since the rescaling will take the UV coupling to zero at 
finite p. Hence we only have one free parameter to work with. Thus we see that it is not 
even guaranteed that that there exists an ai and hi that will connect the UV solution to 
the IR conformal point. In fact, without some symmetry, it seems highly unlikely. Even if 
it is possible to set the ai to zero, it is at best a horrendous fine tuning problem. 

Let us thus accept that the [/F solution attaches to an IR solution with the asymptotic 
behavior in ( p.l3|) , with nonzero ct^. In section 5 we will investigate the consequences of 



this. Before doing this let us round out this section by discussing the most general case 
where g'{To) ^ 0. In this case T will flow as p is varied and thus we have 4 coupled 
nonlinear differential equations. If g'{To) < then T flows toward zero as p is increased. 
One possible IR solution was discussed in |3^] where ^ and 77 increase logarithmically, (p 



increases as a log of a log, and T relaxes to zero. The nice thing about this scenario is that 
the behavior of g{T) and /(T) are known as T ^ 0. Unlike the previous case, there appear 
to be enough integration constants in the UV to match solutions. The generic solution is 



still of the form (|3.13|) with the additional equation 

T ^ asp +To. (3.15) 
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The constraint equation restricts the a^ as 

1 o 1 



2«o + 2«i + 4«3 - 5a^ = (3.16) 

The leading correction to T is yfp(T^p^+ + T-p^-)^ where 

/3± = ^(1 ± v/l-4r(To)//(To)). (3.17) 

If f"(To) > /(To)/4, which is the case for the function in ( p.4| ), then we can write down 
an oscillatory solution with an amplitude and a phase shift. Hence, naively anyway, we 
appear to have enough integration constants to match solutions. But matching to this 
particular IR solution still requires a tremendous amount of fine tuning. 

If /"(To) < /(To)/4 then (3± would both be positive real and there would still be two 
integration constants to adjust. This behavior would also fix the problem of instabilities 
due to the negative tachyon mass near the UV fixed point. At this point, the space is 
locally AdS^ and so the condition for stability is that the tachyon mass satisfy m^ > —4 
P^ , P^ ,|3[] . This is precisely the condition that /"(Tq) < /(To)/4. If this condition were 
satisfied then at the UV fixed point the tachyon would couple to a relevant operator. 

More generic is the case where T relaxes to zero or some other constant Tm while 
the other variables have the asymptotic behavior in ( |3.13| ) . In this case, T — Tjfi falls off 



exponentially as p — * oo. We can see from the equations of motion for T that the second 
derivative starts out positive, but becomes negative at some point before T = 0. Hence, 
it should be possible to find a solution where T starts upward but then relaxes to some 
finite point, while all the other fields grow linearly in p. 

Finally, let us consider g'{To) > 0. Now T fiows away from T = and there appear to 
be no possible IR solutions except for those of (|3.13|) with ct^ > and ( |3.15|) with a^ < 0. 



4. Possible Validity for a Class of Infrared Solutions 



In order to trust the solutions with the asymptotic behavior in (|3.13|) and (|3.15| ), it 



is necessary to check that the higher a' corrections are small. At first glance, this would 
appear to be a miserable failure. The curvature in the Einstein metric ds^ = e~^''^ds'^j 
where ds^ is the metric in (|2.1[), approaches 



R ~ /1^^2_5^^2\ (5a,-ai/2V_ ^^-^^ 
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This blows up as p — > oo because of the constraint in (|3.16| ) . 



The leading string corrections are of the form a' e ^'f'/'^W, where VF is a combination 
of four contracted Riemann tensors. Hence the naive dependence for the a' terms using 
the solutions in (|3.13| ) is 



Using the conditions in ( p.l6| ), one immediately sees that this term grows faster than R, 
and hence the classical approximation breaks down. 

However, W has a field redefinition ambiguity and can be written strictly in terms 
of the Weyl tensor C^jj\s. For the Einstein metric, the nonzero components of the Weyl 
tensor have the form 

C^.^. = g^^^^g.Mi^ + A2g)e'>^-^'\ (4.3) 

where Ai and A2 are constants (different components have different constants), and JF and 
Q are defined as 

jr ^ 2(^-77)77 + I - ij, ^ = ^2 _ ^2 ^ I _ ^ _ 4g-4^_ ^44^ 



It is clear that for a general solution in ( |3.13[ ) , both JF and Q are of order unity, and 
so e~^''^'^W blows up as in (|4.2|) . However, if 02 = ai, then both terms are suppressed. 
We actually should have anticipated this, since now the Einstein metric approaches that 
of the product space Rj^ x M4 x 5*5, up to a conformal transformation. Recall that the 
metric for the 3(4) dimensional Witten model approaches the metric for the product space 
-R5 ^ S5 [Rq X 5*4) near the black hole horizon. 

However, simply setting a2 = cti is not sufficient to insure that e~^'^''^W blows up 
slower than R. From the equations of motion in ( |2.6| ), we see that ^ has corrections of 
order e~2«ip and tj has corrections of order e("o+"i~^*^"2)p/2_ Lg^^jj^g q^^ — q^^^ ^g gj^j jj^ 
general that JF ~ e-2«ip ^nd Q ~ e~2°ip^ and thus 



e 



-3^/^W ~ e(10ai-3ao/2)p_ ^^45^) 



/3 



While the behavior has been softened, it is still not enough to prevent the a' correction 



from dominating i?, since the maximum value for ckq consistent with (|3.16| ) is cto = 3ai. 
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However, there are solutions where W grows even slower than (^4.5| ). If T approaches 
a constant Tm as p -^ oo, then as = and the asymptotic solution has the relations 

cto = ScKi, 0:2 = Oil- (4-6) 

For this case, ^ has corrections of order e~^"i'' and e~^°'^^ while ry has corrections of order 
g-3aip It is a straightforward excercise to check that these corrections cancel in both T 
and Q and therefore 

Hence, not only does the term grow slower than i?, it actually falls off exponentially! 

We can also check the behavior for other types of terms. One class of terms has the 



form [^1 



falling off even faster than the first string correction. There are also terms of the form 

which again falls off rapidly. Terms involving derivatives of the tachyon field also have 
exponential suppression since T approaches a constant value with exponential falloff. The 
behavior of all these corrections leads us to speculate that there is a background that is 
an exact solution of the a model which has the asymptotic behavior defined by (|3.13D and 



There is some fine tuning involved in connecting the desired \JV solution to the IR 
solutions in ([4.6|) . We need to tune such that a\ = 02 and 03 = 0. If g'{To) = then 
CKs = is automatically satisfied. In either case there are enough available integration 
constants on the UV side which can be adjusted so that ( |4.(j| ) is satisfied. Hence, we 
believe that there is a f/F solution that connects to the IR solutions in ( p.l3 ) and (4.6). 



5. Confinement in the Infrared 



In this section, we argue that the IR solutions of ( |3.13|) and (|3.15|) lead to a linear quark 



potential, magnetic screening and a discrete spectrum. The linear quark potential requires 
the additional condition cti > ao- Happily, the solution in ( |4.6| ) satisfies this condition. It 
is possible that other solutions besides (|4.6|) are valid, perhaps there is enough ambiguity 
in the effective action to make their a' corrections small as well. To allow for this, we will 
discuss the general case in this section. 
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5.1. The Quark Potential 



The quark potential is computed using the Wilson loop calculation of [^|36[| . The 
Nambu-Goto action for a string in a curved background is given by 

S = ^ f daSr^detiGMNdc^X^dfsX^). (5.1) 

The quark potential is computed from the rectangular Wilson loop, with two sides along 
the time direction and the other two along one of the spatial directions. Using the metric 
in ( p.l| ), we find that the energy for the quark-antiquark pair is 

(5.2) 



(5.3) 



Chang 


ing variables such that du — —e^'^~^'^>' 
f{u) has the asymptotic behavior 


\dx/ 
^dp, we have 


where 




/(")/«. 



^4 



/(^) ~ 71 ^2 " » 1 

(lOgMj^ 

f[U) ~ u' 1=^ U«l. 

5^2 — CXq 



(5.4) 



Assuming that all ai > 0, the constraint equation ( ^.161 ) implies that the exponent for the 
small u behavior satisfies 7 < 2. li ai > ao, then < 7 < 2, and thus for large enough 
X, the minimum energy configuration consists of a string starting at m = 00 and going 
straight down to the origin at w = 0, traveling a distance x along w = 0, and then going 
back out to w = 00. Since f{u) = at w = 0, there is no cost in energy to separate the 



quarks further, thus this would correspond to electric screening [PO| , p[ . 

If instead ai = ao, then u'' is no longer zero at the origin. In this case, for large 
enough x the string will approach the origin and the quark potential becomes linear. If 
«! < ao, then u'^ blows up at the origin. In this case there will be some nonzero uq such 
that f{uo) is a minimum. For large enough x the string will approach this minimum value. 
Since f{u) is clearly positive for nonzero u, this too will lead to a linear quark potential. 
Assuming that the ai are all of order 1, then the string tension will be of order l/\/Q in 
the units used here. 

If cti < cto, then the coupling seen by the string is in some sense bounded in the 
infrared. For large quark separation, the minimum energy configuration has u > uq, thus 
the coupling approaches e'^o+"i"". This is not to say that the behavior of the metric and 
coupling for u < uq will not affect other physical results, as we will see in the following 
subsections. 
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5.2. Magnetic Screening 

If there is confinement, then there should be a corresponding screening of magnetic 



charge. The computation is along the lines of [0. In particular, we wish to compute the 
potential between a monopole anti-monopole pair. This is accomplished by calculating the 
rectangular Wilson loop for a D string. In type theory, since there are two types of R-R 
Kalb-Ramond fields, there are also two types of D strings. However, only one type of D 
string can end on the electric D3 brane. It is this string that describes the Wilson loop for 
a monopole. 

The calculation is almost identical to the Wilson loop calculation for the fundamental 
string. The only difference is that the integrand in (|5.1| ) is multiplied by a factor of e~'^. 
Hence, the potential is given by 



1 /• . / , ^jdpy 



This time changing variables such that du = e^~^~^'^'/'^dp^ we have 



where now f{u) has the asymptotic behavior 

f{u) ~ u^{\oguY u»l 

^ 2(ai+c^o) ^^ ^ (5-7) 

5^2 + ao 

Since all ai > and 5a2 > ai, then < 7 < 2. Hence, the monopole and anti-monopole 
will screen at large distances 

5.3. WKB Estimates for the Gluehall Spectrum 

The correlator (Tr(F^(x))Tr(F^(0))) is related to scalar Green's functions in the bulk 
0,^ . A discrete spectrum for the bulk particles would correspond to a discrete spectrum for 



the scalar glueballs |]T5|. Since we are dealing with string tensions of order 1, the green's 
functions of massive string states could be important. More to the point, the tachyon 
couples to Tr(F^) [p2| , |l^ . If the tachyon is a bad tachyon, with plane wave normalizable 



states in the bulk, then it cannot have a discrete spectrum. Even if it is well behaved with 
/"(To) < /(To)/4, one might find that the lowest lying state in the spectrum is tachyonic. 
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The linear equations one derives for the smaU fluctuations are hopelessly entangled 
through the tachyon mass term. As a first approximation, one can ignore the cross terms 
and simply consider the diagonal equations. Let us make the usual ansatz and assume 
that the 0"'""'" glueball spectrum is determined by the spectrum for small fluctuations of 
the dilaton equation of motion. Consider dilaton solutions of the form cj) = i/j{p)e^^'^, 
where /c^ = — M^. Then the dilaton equation of motion reduces to 

^v^e-2'^^^''— V; + M^^e-^^g-^ij + -g{T)^i^ = 0, (5.8) 

where cf) in ( ^.8|) refers to the background dilaton fleld. Substituting the metric in (^3) 
into (^]^) and using the equation of motion for the background dilaton fleld, one flnds 

8'^ 1 

tP + M^e^-^^tP (j)"tP = 0. (5.9) 



Deflning u as e^~^^ = (5w^a/2//(To)/32, the equation of motion becomes 

— F{u)—ij + M^QH{u)i; + G{u)ij = 0, (5.10) 

ou ou 

where F{u), H{u) and G{u) have the asymptotic behavior 



(5.11) 



Note that the substitution for u requires that 5a2 > ^i in order that u — *> as p — i> oo. 
Luckily, this inequality is guaranteed by (|3.16D . 



F{u) ^ u^ 


H(u) ~ u 


Gin) . 


'±U 


w >> 1 


logw 


F{u) ~ u 


H{u) ~ u^ 


G{u) ~ w^/2 


u « 1. 



Using the arguments in [^, one can argue that the spectrum is discrete. In fact, we 
can say more. Since we know the asymptotic behavior, one can do the WKB calculation 



described in [^5|. The asymptotic behavior is similar to that of the three dimensional 
Witten model [jl5|, although there will be log corrections because of the UV behavior 
of G{u). The G{u) term does not affect the IR turning point to next to leading order. 
Borrowing the results in [|25| and inserting the G{u) correction, we flnd that the WKB 0"'""'" 
glueball spectrum is 

„. ^ ^cywrn „/^^, 2 \^/^^\ 

Q V log[Cm(m + l)]y \log^mJ 

(5.12) 

18 



where C is a constant of order unity. Corrections to (|5.12|) coming from the inclusion of 
cross terms are of order mj log m. The log factors in ( p.l2|) have the effect of dragging 
down the glueball masses, which is clearly caused by the running of the coupling. This is 
especially true for the lightest state. If C = 1 and one assumed that ( |5.12| ) were exact, then 
one would find that the lightest state is actually tachyonic. At this point we are unable to 
prove that a tachyonic state does not appear. We can say that if a tachyonic state were 
to exist for the dilaton equation in ( |5.8|) , then a minimum requirement is that there is 
some region where p^(\)" < —1/2. One can see that this is never true if we only consider 
the leading order term in ( |2.15| ) for (p. Numerical analysis seems to show that it is not 
true if the higher order terms are also included. Once we are safely in the infrared, then 
(p" will be exponentially suppressed and so p'^cp" > —1/2 in this region as well. However, 
inclusion of the cross terms might drop the masses enough to leave a tachyonic glueball 
state. Hopefully this does not occur. 

Since Q ~ A^, we see that the glueball masses are a factor of A^/^ smaller than 
the square root of the string tension. This is not as severe as the 4 dimensional Witten 
model, where the masses are a factor of A^'^ smaller |^^. We can also compare the 



ratio of the first excited glueball to the lowest mass state. Since our knowledge of C is 
limited, we can only put a lower bound on the mass ratios, which is reached when C —>■ oo. 
Plugging in numbers, we find that this ratio satisfies mi^_,_/m_|_+ > -\/3 ~ 1.73. The WKB 
approximation for this ratio in the 4 dimensional Witten model is a/8/3 ~ 1.63 |2^. The 
SU(3) lattice result is 1.77 ± .14 JH 



We can also make WKB approximations for the glueball states with nontrivial 5*0(6) 
quantum numbers j^^. If we write the dilaton field as (p = ip{p)e^^'^Yi{Q^), then the 



dilaton equation of motion reduces to 



d' 



tP + M^e^-^'^tp -l{l + 4)e-^'^tp--(p"ip = 0. (5.13) 



Qp2^ ' ^ V ■ ^ ^2 



If cti > a2, which is the case for ( [4.6| ), then the M^ term will dominate the i{£ + 4) term as 
M — i> 0. In this case, the WKB calculation is a straightforward generalization of the / = 
result and is 

_2 4Ca/2/(To) / , 2 \ .. f rn \ 

Q V log[Cm(m + l)]y Viog^^^y 

(5.14) 

If ai < a2, then the asymptotic behavior changes near u = 0, but this effect only leads to 

corrections of order 1 in m. 
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We complete this subsection by computing the WKB spectrum for the '" gluebaUs. 
In the Witten model, this was done by studying the wave function for the 5 dimensional 
dual of the NS-NS two form field, which is a vector [1T9|J2^ . The component along the 
Euclidean time direction is a scalar in the four dimensional theory and is odd under parity. 
In the case that we are studying, there is no euclidean circle, so we need to consider instead 
the wave function for a scalar field. The scalar field is the IIB axion which is an R-R scalar. 
Actually, there are two scalars, but only one couples to Tt{FF) on the electric D3 branes. 
The other scalar couples to Tt{FF) on the magnetic D3 branes. 

Since the axion comes from the R-R sector, it does not couple to the dilatona. Hence, 
the wave function has the form 

op op 



Defining a{p) = e '^'(p{p), ( p.l5| ) can be rewritten as 

92 



dp' 



t/j + M^e^-^^il; - {(f)" + (f)'^)ij = 0. (5.16) 



As was the case for the dilaton, the (p terms do affect the next to leading order contribution 
of the t/F turning point. The (J)' term also changes the next to leading order contribution 



of the IR turning point. Using the analysis in [gSj], one finds that the WKB mass expression 
for the '" states are 



M — ^^-— m { m + 1 + — 7— — ; —T + O — 7j — m > 1 



Q V 5^2 - «! log[Cm(m-|- 1)]/ Vlog^m 

(5.17) 

where C is the same constant as in ( |5.12|) . Unlike the 0"'""'" states, but like the '" states in 



the Russo generalization of the Witten model p3| , p4| , |27|| , there is model dependence in these 
glueball masses. Using the solutions in ( [4.6| ) we find as a lower bound for the mass ratio 

m ^/m-|-+ > v7/2 ^ 1.32. The most recent lattice results have m ^/m++ = 1.78 ± .24 

for SU{2) and m-+/m++ = 1.50 ± .04 for SU{3) gH 



The axion does couple to the tachyon [28|, but this does not affect the WKB approximation 
to next to leading order. 
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6. Discussion 

One of the main differences between tlie results found iiere and results found for tlie 
Witten model is that the type model has a QCD scale set by dimensional transmutation. 
Hence, the glueball mass is directly related to the string scale. We have also found dif- 
ferences for the mass ratios. The principle reasons for these differences is that the type 
model is four, not five dimensional in the UV, and because the type model has a running 
coupling constant. The A^ dependence in the mass to string scale ratio is also affected by 
the differences in dimension. The Witten model has a factor of A^~^/^, while the type 
model has a factor of N~^''^ 

The existence of the confining solutions does not contradict other results for IR fixed 
points. The fine tuning involved to reach such solutions could correspond to tuning mass 
parameters on the field theory side, so that for instance, the six scalars are massless after 
including their quantum corrections. In this case, the field theory could have the IR fixed 
point described in ||3^. The confining solutions, which do not need to be fine tuned, have 



massive scalars, and so are in the same universality class as pure SU{N)B 

We have seen that while confinement in the infrared is generic, there exists solutions 
where many of the string corrections vanish exponentially. One such solution has a con- 
stant tachyon expectation value. In some sense, this is closer to the spirit of Polyakov's 
noncritical string for QCD, where the nonzero tachyon expectation value, "peacefully con- 
densing in the bulk" , plays the role of a nonzero central charge [^ , ^ . On the other hand. 



the five form flux through 5*5 clearly plays an important part here; without it we would 
not have the log scaling in the UV. 
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